An exactly solvable model of a highly efficient thermoelectric engine 



Martin HorvalE] and Tomaz Proser{l] 
Faculty of mathematics and physics, Department of Physics, 
University of Ljubljana, Jadranska 19, SI- 1000 Ljubljana, Slovenia 

Giulio CasatH] 

Center for Nonlinear and Complex Systems, Via Vallegio, 11, 22100 Como Italia 

(Dated: June 18, 2009) 

We propose a simple classical dynamical model of a thermoelectric (or thermochemical) heat 
engine based on a pair of ideal gas containers connected by two unequal scattering channels. The 
model is solved analytically and it is shown that a suitable combination of parameters can be chosen 
such that the engine operates at Carnot's efficiency. 
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In the frame of non-equilibrium thermodynamics an 
heat engine is a machine generating work while exchang- 
ing heat with two heat baths at different temperatures 
Ti and T2. The usual goal in a construction of realistic 
heat engines is to increase the efficiency as far as possible 
towards the theoretical upper limit 77carnot = 1 — ^1/^2, 
assuming T2 > Ti. We are here interested in an en- 
gine without moving mechanical parts, i.e. which could 
operate in a non-equilibrium steady state^ such as for ex- 
ample a thermoelectric or thermochemical couple. Such 
an engine - or a refrigerator if the operation is reversed - 
would have immense practical advantages over piston or 
compressor based engines for obvious reasons, also due 
to possibilities of drastic miniaturization [1 . 

Here we present an abstract model of a heat engine 
that can mimic the essential features of a realistic heat 
engine based on the thermoelectric effect [11 [3l [4] and 
which can be treated and solved analytically. It is based 
purely on deterministic classical dynamics and stochas- 
tic baths. The model is composed of two thermochemical 
reservoirs of ideal gas of equal point particles connected 
by two one-dimensional wires indexed by z G {1,2}. In 
the middle of each wire we place a deterministic and en- 
ergy conserving scatterer, which either reflect or transmit 
the particle depending on its kinetic energy e. This be- 
haviour is completely described by the transmission func- 
tion ri{e) G {0, 1} of the i— th scatterer. We use units in 
which particle mass m, particle charge e and Boltzmann 
constant /cb equal m = e = = 1. In this paper we 
show that in the steady state a non-vanishing circular 
particle current exists only if the transmission functions 
are energy dependent. Then we show that for a suit- 
able combination of parameters the engine operates in a 
reversible way with the Carnot's efficiency. 

The scheme of the heat engine is shown in figure [l] In 
the wires we introduce bias forces Ei (say electric fields), 
which can be described by bias voltages Ui or any other 
form of external potential energy which can be used to 
extract useful work. In the stationary state, at some tem- 
perature difference, there is a non-zero (circular) particle 



current in the wires that, by climbing against the elec- 
tric potential, can perform useful work. In the left (right) 
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FIG. 1: Schematic figure of the heat engine. The possible de- 
terministic scattering mechanisms are depicted symbolically. 

reservoir, the particles are at chemical potentials /iL (mr) 
and temperature Tl (Tr) (here we assume Tl > Tr) and 
are effused into the wires with the injection rates pi7L 
iPilK) into the first i = 1 and the second i = 2 channel, 
respectively, where Pi G [0, 1], pi + P2 = 1, represent the 
relatives openings into the two channels. We note that 
both in the reservoirs and in the channels the motion of 
particles is assumed to be (quasi) one dimensional^ so 
we consider a single component of the velocity. The in- 
jection rates are connected to the chemical potentials jiy 
and inverse temperatures jSi, = l/Tj^ via the formula 



liypy = log(C/3^7^) , G {L, R} , 



(1) 



where C is a constant depending only on properties of 
particles and on geometry of the reservoir opening [4]. 
The velocity v of effused particles is distributed in each 
side according to a canonical distribution 



(2) 



where 6{v) = {1 : v > 0;0 : otherwise) is the unit step 
function, and (Jl = l,crR = — 1. In the steady state, 
the particle currents jp^i within i-th wire and the heat 
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currents [2^ exchanged with the i^-side bath and i-th wire 
jq^ilj, are given by 

jp,i = Pi(7L^L,i - 7R^R,i): (3) 
jq,ik = PiilhQL^i - 7R(te,i + ^R,i^i)), (4) 
jq,i\R = Pi(7L(<7L,i - ti^,iUi) - 7R<?R,i). (5) 

Here we have introduced the transmission probabihty tj^^i 
for a particle to transit from the i^-side to the other side 
over the i-th wire, and its average kinetic energy ^j.^^, 
exphcitly defined in terms of the first two statistical mo- 
ments of the energy distribution of the effused particles 
transmitted through the zth wire 

e-'^'n(e-^) (l,e), (6) 
+ y) (7) 



(^L,i,g'L,i) = de 

J max{0,t 

(tR,i,gR,i) = pR dee"^^Vi(e 

Jmax{0,-Ui} \ 



The terms Ui/2 in the arguments of transmission func- 
tions Ti imply the assumption of linear potential and scat- 
terers being in the middle of each wire. However different 
assumptions (say of bias potential steps at the left/right 
of each scatterer) could be treated straightforwardly. By 
imposing the condition of stationarity jp^i + jp,2 = 0, we 
obtain, from eq. (3), the injection rates 7l and 7r 

^ 7L ^ 7R ^ Pljth,! - fa,l) +P2ftL,2 - fa,2) .g^ 
7 7 Pl(^L,l +^R,l) +P2(^L,2 +^R,2) ' 

where 7 = (7L +7r)/2 [5j. The resulting particle current 
in the wires 



jp,i = 27P1P2 



^L,l^R,2 — ^L,2^R,1 



Pl(^L,l +^R,l) +P2(^L,2 +^R,2) ' 



(9) 



determines the working power P = jp^iiUi — U2) while 
the ingoing heat flux is equal to 



Q = jg,l|L + jg,2|L 



Pl(^L,l + ^R,l) + P2(^L,2 + ^R,2) 



(10) 



where we have introduced the ratio Su,i = qu,i/tu,i (which 
is shown below to be connected to the Seebeck coef- 
ficient). The efficiency of the heat engine is then de- 
fined as 77 = P/\Q\. Notice that the particle current 
jp^i, and so also the power P, is proportional to the 
determinant of the matrix of transmission coefficients 
D = det{U^i} = tL,i^R,2 - ^L,2^R,i. The optimal per- 
formance of the heat engine for a given configuration 
of temperatures and scatterers, is obtained by finding 
appropriate fields Ui and U2 which maximize r]{Ui,U2). 
This can be done, in general, only numerically since it 
requires a solution of coupled transcendental equations. 

Note the following important observation: D = 0, and 
hence the currents vanish, despite non- vanishing temper- 
ature difference, if the scatterers are energy independent 



ri{e) = const. This fact is a simple consequence of time- 
reversal properties of individual deterministic trajectories 
which connect the two baths and remains valid for scat- 
tering channels in higher dimension (e.g. like in Ref. [4|). 

Our ideas are demonstrated in a heat engine with bias 
voltage only in the first wire {U2 = 0) and for the simplest 
nontrivial, step-like, transmission functions 



ri{e)=e{s,{e-ei)). 



(11) 



The direction of the steps, at the energy thresholds e^, is 
determined by the signs Si G {1,-1}. Simple mechanical 
realizations for both signs Si G {1,-1} are schematically 
depicted in fig. [T] We numerically determine the poten- 
tial Ui that maximizes the efficiency for a given configu- 
ration of scatterers. The optimal efficiency 77* = max^/^ rj 
and the corresponding power P* are shown in figure [2j 
In the case Si = S2 = —1 the scatterers transmit only 
slow enough particles. The regions of high power and 
high efficiency overlap and are positioned almost sym- 
metrically near the axes. The exact symmetry is broken 
because the electric field is only applied to the first wire. 
The cases si = —S2 = —1 and si = —S2 = 1 describe 
a similar situation, where one scatterer transmit the fast 
particles and the other scatterer transmit slower ones. 
This case is the most efficient and rj here may nearly ap- 
proach r^carnot- Howcvcr, as expected, the regions of high 
efficiency and high power only slightly overlap. In the last 
case Si = 52 = 1 the scatterers only transmit fast enough 
particles. The region of high efficiency is located parallel 
to the line ei = 62- The highest power is obtained for 
energy steps matching the baths temperatures ei = Tl, 
62 = Tr. a detailed analysis of the relaxation process 
shows that the convergence time to the non-equilibrium 
steady state is strictly finite for non- vanishing bias po- 
tentials, and is given by r* ~ 1/ min{|/7i|, \ U2\}. 

In the following we show that the results drastically 
simplify in the linear regime of small relative temperature 
difference. In this regime the meaningful bias potentials 
are also small and we may approximate the exact particle 
and heat fluxes in the wires with their linear expansions 
in the temperature difference Tl — Tr, injection rate dif- 
ference (^7 = 7r — 7l, and potentials Ui. Expressing Sj 
with the chemical potential difference (^/i = /iR — /iL, set- 
ting /iR + /iL =0, we write the particle and heat fluxes 
in the linear response limit as 



Jp,i 
Jq,i 



--Pil{-p9i{Ui 
■- Pij{-phi{Ui 



5^) - 



(12) 
(13) 



where P = {Pl ^ /5r)/2 and 5p = Pr - Pl- The ex- 
pansion coefficients gi^ hi and ki are statistical moments 
of a canonical energy distribution of particles that are 
transmitted over the i-th wire 



(gi^hi.ki) = p 
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det 



-fte 



n(e)(l,e,e2) (14) 
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FIG. 2: The relative power P*/'/7carnot (a) and the relative 
efficiency ?7*/?7carnot (b) as functions of the energy thresholds 
Ci for different sign configurations (si,S2) (indicated on the 
left) at unit mean injection rate 7=1, equal channel openings 
Pi = 1/2, and bath temperatures Tl = 2 and Tr = 0.5. 



and depend only on the transmission function and tem- 
perature. Notice that in the linear response limit the 
heat fluxes at the left and right side are equal in contrast 
to the general (non-linear) case Q, ([5|. The coefficients 
gi and hi represent the average transmission probability 
of particles across the i-th wire and their average en- 
ergy at zero bias fields. Instead of hi and ki it is more 
convenient to work with the average energy per particle 
Si = hi/gi and the coefficient Ki = ki — giSf. Note 
that Ppijgi, P'^pijKi and pSi can be interpreted as the 
particle conductance^ the heat conductance and the See- 
beck coefficient^ respectively. By imposing the stationar- 
ity condition jp^i + jp^2 = we obtain the difference of 



the chemical potentials between baths: 

^{pigiSi +P2^2*5'2) ^PigiUi +^2^2^/2 



5fi = — - 



Pigi +P2^2 



(15) 



where ^ = —SP/P is the relative temperature difference 
which is related to the Carnot efficiency r^carnot = 1^1- 

Let us now introduce the auxiliary quantities: dif- 
ference of energies per particle in the two wires AS* = 
S2 — Si^ difference of potentials SU = U2 — Ui^ trans- 
mission probability through both wires G = {{pigi)~^ -\- 
(^2^2)"^)"^, and the figure of merit of the heat engine 
efficiency [6^ 



PiKi + P2K2 



(16) 



We can now write the particle current in the first wire 
jp^i and the ingoing heat flux Q = jq^i + jg,2 elegantly as 



Q 



(B-yGiASSU 



pjG{-5U - 



■^AS), 

■ 1/y)) , 



(17) 
(18) 



whereby the power and the efficiency of the heat engine 
are P = jp^iSU and 77 = P/|Q|, respectively. Notice that 
all expressions just depend on the potential difference 
6U. In the linear response regime the potential that max- 
imises the efficiency can be found analytically by solving 
the equation dr]/d{5U) = 0. The explicit solutions are 



(5C/* = CA5(l-(v/rT^-l)/t/), (19) 
= e7/3GA5( vTTi; - l)/y , (20) 
Q* = -ij(3G{ASf./lT^/y , (21) 

yielding the optimal efficiency (equivalent to Eq. (14) of 
L6) 



V 



^carnot 



(i + H(i-ym;)) 



(22) 



with the corresponding power P* = jp^iSU^. The rela- 
tive optimal efficiency f]* /rjcamot ? as expected, depends 
only on y and is monotonic in the latter. Therefore it is 
meaningful to treat y as the figure of merit of heat engine 
efficiency. We note that if the transmission functions are 
given by ( 11 ), the results ( T9p2 ) are explicit as all the ex- 
pressions are explicit rational functions of the moments 
of the Laplace transform of the transmission functions 



(14), which in turn are simple algebraic functions of 
and exp(— /3e^). From the equations (19) and (20) we 
can recognize that P* oc ^carnot ^ consequently, in the 
linear response regime the power-output is rather small. 
The optimal efficiency and the corresponding power as 
function of ci and €2 in linear response regime are shown 
in figure [3] They are quite similar to those obtained in 
the non-linear regime shown in figure [2] The important 
difference between nonlinear and linear regime results is 
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FIG. 3: The density plot of the relative power 10P*/?7, 
(a) and the relative efficiency ?7*/?7carnot (b) in the linear 
response regime, at the bath temperatures Tl = 1.01 and 
Tr = 0.99. (Other details the same as in figure [5]) 



si = —52 = 1 (and similarly for si = —52 = —1) we 
can reach the Carnot efficiency in the limit ri ^ oo 
following the curve x 2exp(— ri) along which the 
efficiency algebraically increases as 77* x 77carnot(l — 
2r]"^), and the power exponentially decreases as P* x 

In conclusions, we have proposed a simple exactly solv- 
able classical-mechanical model of thermoelectric (or bet- 
ter to say, thermochemical) heat engine. We presented 
closed form solutions for the steady state of the engine 
in linear and non-linear regimes. A variable thermody- 
namic efficiency has been found, as a function of the sys- 
tem's parameters, which can become arbitrary close to 
Carnot's in an appropriate regime. 

Finally we would like to draw the reader attention to 
the following point: it is possible to argue that our model 
is quite abstract in nature and therefore far from possi- 
ble realistic implementations. We think on the contrary 
that this is the main advantage of our approach. After 
more than 50 years during which thermoelectric efficiency 
did not substantially increase we propose here to take a 
completely opposite point of view. Starting from funda- 
mental microscopic equations and considering the most 
general schematized framework, we hope to understand 
the basic dynamical mechanisms which can lead to an 
increase of thermoelectric efficiency. In this spirit, the 
model discussed here is a step in this direction. In ad- 
dition our model should be relevant for a theoretical de- 
scription of nanoscopic heat engines, for example a pair of 
thermoelectrically coupled quantum dots. However, our 
model would be a good approximation to the real system 
only in a rather restricted situation of (i) non-interacting 
charge carriers, (ii) negligible phonon contributions to 
heat transport, (iii) coherence length longer than wires, 
which (iv) should be quasi one-dimensional. 



that the latter only depends on the difference dU of bias 
potentials and can be made temperature independent by 
expressing the energy steps in the transmission function 
Ti{e) with the parameters = . Consequently, the 
power P* and the efficiency 7^* as function of (61,62) in 
the cases si = —52 = 1 and si = —52 = —1 are exactly 
symmetric w.r.t. exchange of the parameters e^. These 
cases are also the most efficient. The region of high effi- 
ciency is squeezed towards the 62 or ei axes in the cases 
5i = —52 = 1 and 5i = —52 = — 1, respectively. 

We have performed exact analytical calculations of ex- 
pressions ([19][22| for the case of equal channel openings 
Pi = 1/2. In the cases 5i = 52 = 1 and 5i = 52 = —1 
the maximal efficiency r^max = niax^^^^2>o ^* is reached 
at finite (ri,r2) and is equal to 77max = 0.066 T^camot and 
^max = 0.091 T^carnot, respectively. However, in the case 
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